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Abstract 

We investigate the Heston model with stochastic volatility and exponential tails as 
a model for the typical price fluctuations of the Brazilian Sao Paulo Stock Exchange 
Index (IBOVESPA). Raw prices are first corrected for inflation and a period span- 
ning 15 years characterized by memoryless returns is chosen for the analysis. Model 
parameters are estimated by observing volatility scaling and correlation properties. 
We show that the Heston model with at least two time scales for the volatility 
mean reverting dynamics satisfactorily describes price fluctuations ranging from 
time scales larger than 20 minutes to 160 days. At time scales shorter than 20 min- 
utes we observe autocorrelated returns and power law tails incompatible with the 
Heston model. Despite major regulatory changes, hyperinflation and currency crises 
experienced by the Brazilian market in the period studied, the general success of 
the description provided may be regarded as an evidence for a general underlying 
dynamics of price fluctuations at intermediate mesoeconomic time scales well ap- 
proximated by the Heston model. We also notice that the connection between the 
Heston model and Ehrenfest urn models could be exploited for bringing new insights 
into the microeconomic market mechanics. 
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1 Introduction 



In the last decades the quantitative finance community has devoted much 
attention to the modeling of asset returns having as a major drive the im- 
provement of pricing techniques [1] by employing stochastic volatility models 
ameanable to analytical treatment such as Hull- White [2], Stein-Stein [3] and 
Heston [4] models. 

Despite differences in methods and emphasis, the cross fecundation between 
Economics and Physics, which dates back to the early nineteenth century (see 
[5] and [6]), has intensified recently [7]. Following the statistical physics ap- 
proach, substantial effort has been made to find microeconomic models capable 
of reproducing a number of recurrent features of financial time series such as: 
returns aggregation (probability distributions at any time scale) [8,9], volatil- 
ity clustering [11], leverage effect (correlation between returns and volatilities) 
[12,13], conditional correlations [14,15] and fat tails at very short time scales 
[7]. A central feature of economical phenomena is the prevalence of intertwined 
dynamics at several time scales. In very general terms, one could divide the 
market dynamics onto three broad classes: the microeconomic dynamics at the 
time scales of books of orders and price formation, the mesoeconomic dynam- 
ics at the scales of oscillations in formed prices due to local supply and demand 
and the macroeconomic dynamics at the scales of aggregated economy trends. 

The literature on empirical finance have emphasized the multifractal scaling 
[16] and the power law tails of price fluctuations. However it has been shown 
[17] that very large data sets are required in order to distinguish between a 
multifractal and power law tailed process and a stochastic volatility model. 
In this paper we, therefore, deal with the mesoeconomic dynamics as it would 
be described by the Heston model with stochastic volatility and exponential 
tails. 

Recently, a semi-analytical solution for the Fokker-Planck equation describing 
the distribution of log- returns in the Heston model has been proposed [8] . The 
authors were able to show a satisfactory agreement between return distribu- 
tions of a number of developed market stock indices and the model for time 
scales spanning a wide interval ranging from 1 to 100 days. More recently, 
the same model has also been employed to describe single stocks intraday 
fluctuations with relative success [10]. 

In this paper we show evidence that the Heston model is also capable of 
describing the fluctuation dynamics of an emerging market price index, the 
Brazilian Sao Paulo Stock Exchange Index (IBOVESPA). We employ for the 
analysis 37 years of data since IBOVESPA inception in January, 1968 and 
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approximately four years of intraday data as well. In this period the Brazilian 
economy has experienced periods of political and economical unrest, of hyper- 
inflation, of currency crises and of major regulatory changes. These distinctive 
characteristics make the IBOVESPA an interesting ground for general mod- 
elling and data analysis and for testing the limits of description of the Heston 
model. 

This paper is organized as follows. The next section surveys the Heston model, 
its semi-analytical solution and its connection to Ehrenfest urn models. Section 
3 discusses the pre-processing necessary for isolating the fluctuations to be 
described by the Heston model from exogenous effects. Section 4 describes the 
data analysis at low (from daily fluctuations) and high (intraday fluctuations) 
frequencies. Conclusions and further directions are presented in Section 5. 



2 The Heston Model 

2.1 Semi- analytical solution 

The Heston Model describes the dynamics of stock prices S t as a geomet- 
ric Brownian motion with volatility given by a Cox-Ingersoll-Ross (or Feller) 
mean-reverting dynamics. In the Ito differential form the model reads: 

ds t =s tf Hdt + s t ^r t dw (t) (i) 

d Vt = - 1 [ Vt -6]dt + Kyfa dW^t), 

where v t represents the square of the volatility and dWj are Wiener processes 
with: 

=0, 

(dWj{t) dW k (t')) = S(t - t') [5 jk dt + (l- 8 jk )pdt\ (2) 

The volatility reverts towards a macroeconomic long term mean squared volatil- 
ity 9 with relaxation time given by 7 -1 , fit represents a drift at macroeconomic 
scales, the coefficient ^JW t prevents negative volatilities and k regulates the am- 
plitude of volatility fluctuations. 

As we are mainly concerned with price fluctuations, we simplify equation (1) 
by introducing log-returns in a window t as r(t) = \n(S(t)) — ln(5'(0)). Using 
Ito's lemma and changing variables by making x(t) = r(t) — Jq ds /j, s we obtain 
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a detrended version of the return dynamics that reads: 



dx = — ^ dt + y^dWo- 



(3) 



The solution of the Fokker-Planck equation (FP) describing the unconditional 
distribution of log-returns was obtained by Dragulescu and Yakovenko [8] 
yielding: 



+00 



P, M = / ^ e 



dp. 



(4) 



where 



a = 



2 7 # 



. . Qt\ a 2 -r 2 + 2 7 r , /at 

cosh — M ^ smh — 

1 2 J 2 7 fi V 2 



= yT 2 + k 2 (pi - ip x 
T = 7 + ipnp x . 



(5) 

(6) 

(7) 
(8) 



This unconditional probability density has exponentially decaying tails and 
the following asymptotic form for short time t <C 7 -1 : 



„, , 2 1 - a e- x / 2 roT (2ax 2 \ 4 T ^ / 2ax 2 



(9) 



where Kp{x) is the modified Bessel function of order (3. 



2.2 Volatility autocorrelation, volatility distribution and leverage function 



The volatility can be obtained by integrating (1) and is given by: 



1 

Vt = ( Vo - 6) e" 71 * + 6 + kJ dWxiu) 



-7(t-u) 



'V,,.. 



(10) 



A simple calculation gives the stationary autocorrelation function: 



C{t I 7, 9, k) = lim 

t^QO 



(v t v t+T ) - (v t ) (v t+T ) e 



■JT 



ff 1 



a 



(11) 
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The probability density for the volatility can be obtained as the stationary 
solution for a Fokker-Planck equation describing v t and reads: 



U(v) = 



T{a) 6 a ' 



-av/8 



(12) 



implying that a > 1 is required in order to have a vanishing probability density 
as v — > 0. 

The leverage function describes the correlation between returns and volatilities 
and is given by [13]: 



L(r I 7 , 6, k, p) ee lim {dXt ^ } = pn H{r) G(r) e^, 

((dx t+T ) 2 ) 



(13) 



where dx t is given by (3), H{r) is the Heaviside step function and: 



G(r) = 



v t exp 



f J t t+T dW^u) v u 



(vtY 



(14) 



To simplify the numerical calculations we employ throughout this paper the 
zeroth order appoximation G(r) ~ G(0) = O^ 1 . The approximation error 
increases with the time lag r but is not critical to our conclusions. 



2.3 Relation to Ehrenfest Urn Model 



The Ehrenfest Urn (EU) as a model for the market return fluctuations has 
been studied in [18], in this section we observe that Feller processes as (1) can 
be produced by a sum of squared Ornstein-Uhlenbeck processes (OU), and 
that OU processes can be obtained as a large urn limit for the EM. 

To see how those connections take shape we start by supposing that Xj are 
OU processes: 

dX j (t) = ~X j (t)dt + ^dW j (t), (15) 



where dWj describe d independent Wiener processes. In this case the variable 
v(t) = 52j=iXj(t) is described by a Feller process as (1) [19]. The volatility 
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process in (1) emerges from OU processes by applying Ito's Lemma to get: 

d d 2 d 

dv t = -bdtJ^Xj + a X] x i dw 3 + 7E dW j- ( 16 ) 



4 , 



Using the definition of v and the properties of the Wiener processes it follows 
that: 



dv t 



d 2 . 

-a - bv t 



dt + a^W t dW. (17) 



The volatility process in (1) can be recovered with a few variable choices: 
a = k, b = 7 and 9 — A dynamics with vanishing probability of return 
to the origin requires the volatility to be represented by a sum of at least two 
elementary OU processes, equivalently we should have a — | > 1. 

An OU process can be obtained as a limit for an Ehrenfest urn model (EM). 
In an EM N numbered balls are distributed between two urns A and B. At 
each time step a number 1,...,N is chosen at random and the selected ball 
changes urn. We can specify S(n) = (si(n), s N (n)) as the system microstate 
at instant n, where Sj(n) = {±1} indicates that the ball j is inside urn A if 
Sj — 1, we also can define the macrostate M{n) = J2f=i s j( n ), which dynamics 
is described by a Markov chain over the space 0, N with transition matrix 
given by Q(k, fc+1) = 1-k/N, Q(k, fc-1) = k/N and Q(j, k) = if \j-k\ ^ 1. 
An imbalance in the population between the two urns generates a restitution 
force and, consequently, a mean reverting process for M(n). Applying the 
thermodynamic limit N — > oo, rescaling time to t = n/N and introducing a 
rescaled imbalance variable as: 



we recover an OU process: 

dX t = -2X t dt + dW t . (19) 



Using this connection we could speculate on a possible microscopic model 
that would generate a stochastic dynamics as described by the Heston model. 
Supposing that market agents choose at each time step between two regimes 
(urns) that may represent different strategies, such as technical and fundamen- 
tal trading or expectations on future bullish or bearish markets, imbalances 
between populations in each regime would be a source of volatility. In this pic- 
ture, the condition d > 2 would imply that at least two independent sources 
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Fig. 1. A: Inflation discount factor used to adjust past values to January 2005. The 
hyperinflation period (average of 25% per month) from January 1987 to July 1994 
is also indicated. B: Inflation adjusted IBOVESPA index. 

of volatility would be driving the market. We shall develop this connection 
further elsewhere. 



3 Data Pre-processing 



3.1 The data 



Two data sets have been used: IB1 consisting of daily data from IBOVESPA in- 
ception on January, 1968 to January, 2005 and IB2 consisting of high-frequency 
data sampled at 30 second intervals from March 3, 2001 to October 25, 2002 
and from June 6, 2003 to August 18, 2004. 



3.2 Inplits and Inflation 



The dataset IB1 has been adjusted for eleven divisions by 10 (inplits) in- 
troduced in the period for disclosure purposes [20] and also for inflation by 
the General Price Index (IGP) [21]. In figure 1A we show the discount fac- 
tor from date t to current date T, Ir{t), used for correcting past prices S t as 
Sf = Sflrit). Figure IB shows the resulting inflation adjusted prices. The 
hiperinflation (average of 25% per month) period from January 1987 to July 
1994 is also indicated in both figures by dashed lines. 
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Fig. 2. A: Autocorrelation function in the 1968-2005 period showing one day mem- 
ory. B: Historical autocorrelation function measured with 250 days moving windows 
for one day lag. The behavior became compatible with a random walk after the Col- 
lor Plan in 1990. 

3. 3 Detrending 



For our analysis of price fluctuations it would be highly desirable to identify 
macroeconomic trends that may be present in the data. A popular detrend- 
ing technique is the Hodrick-Prescott (HP) filtering [22] which is based on 
decomposing the dynamics into a permanent component, or trend xp(t), and 
into a stochastic stationary component x(t) as r(t) = xp(t) + x(t) by fitting 
a smooth function to the data. Any meaningful detrending procedure has to 
conserve statistical properties that define the fluctuations. However, in our 
experiments we have noticed that the HP filtering may introduce subdiffusive 
behavior at long time scales as an artifact when applied on first differences of a 
random walk. In this paper, in the absence of a reliable detrending procedure, 
we assume that the major long term trend is represented by inflation. 



3.4 Autocorrelation of Returns 



In the period span by data set IB1 the Brazilian economy (see [23] for a brief 
historical account) has experienced a number of regulatory changes with con- 
sequences for price fluctuations. In [23] it has been observed that the Hurst 
exponent for daily IBOVESPA returns shows an abrupt transition from an 
average compatible with long memory (H > 0.5) to a random walk behavior 
(H = 0.5) that coincides with major regulatory changes (known as the Collor 
Plan). We have confirmed the presence of memory by measuring the autocor- 
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Fig. 3. A and C: ranked representations of, respectively, drawdowns and drawups. 
B and D: full lines represent cummulative empirical distributions, dashed lines rep- 
resent minimum and maximum values for 100 shuffled versions of IB3. 



relation function for the entire time series (figure 2A). We have also measured 
the historical autocorrelation for one day lags using a 250 days moving win- 
dow and confirm an abrupt behavior change coinciding with the Collor Plan in 
1990 (figure 2B). As the Heston model assumes uncorrelated returns and this 
feature only became realistic after the Collor Plan, the following analysis is 
restricted to a data set consisting of daily data from January 1990 to January 
2005 (IB3). 



3.5 Extreme Events and Abnormal Runs 



Abnormal runs are sequences of returns of the same sign that are incompatible 
with the hypotheses of uncorrelated daily returns. We follow [24] and calculate 
for the data set IB3 the statistics of persistent price decrease {drawdowns) 
or increase {drawups). We then compare empirical distributions of runs with 
shuffled versions of IB3. In figure 3B we show that a seventeen business days 
drawdown from March 7, 1990 to March 29,1990 is statistically incompatible, 
within a 98% confidence interval, with equally distributed uncorrelated time 
series. Observe that the largest drawup shown in figure 3C correponds to the 
subsequent period from March 30, 1990 to April 5, 1990 and that the Collor 
Plan was launched in March 1990. We therefore expunged from data set IB3 
abnormally correlated runs that took place in March, 1990. 
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4 Data Analysis 



4-1 Low Frequency 



After deflating prices and expunging autocorrelated returns, four independent 
parameters have to be fit to the data: the long term mean volatility 9, the 
relaxation time for mean reversion 7 -1 , the amplitude of volatility fluctuations 
k and the correlation between price and volatility p. It has become apparent in 
[9] that a direct least squares fit to the probability density (4) yields parameters 
that are not uniquely defined. The data analysis adopted in this work consists 
in looking for statistics predicted by the model that can be easily measured 
and compared. In the following subsections we describe these statistics. 



4-1.1 Long term mean volatility 6 

A straightforward calculation yields the second cummulant of the probability 
density (4) as: 



c 2 (t) = (x(tf)-(x(t)Y 



-a 



dk 2 



(20) 



k=0 



■■6t[l-e]+e- l-e" 7 * 
7 L 



where 



K 



7 



P 



Ik 

4 7 



(21) 



As k C 7 we drop terms of O y^J or superior. A non-biased estimate for the 
above cummulant can be calculated easily from data using: 

2 

(22) 



N 



C 2 (t) 



N 



3T§ 



7=1 . 



X 



(i) 



N 



x 



(t) 



where x^ stands for t-days detrended log-returns. The long term mean volatil- 
ity is estimated by a linear regression over the function c 2 (i) as shown in figure 
4A. 
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Fig. 4. A: Linear regression for the second cumulant C2(t) at several time scales. The 
angular coefficient estimates the long-term mean volatility 6. B: Linear regression 
for estimating the shorter relaxation time 7] -1 - The t statistics is shown below the 
estimate and assumes 18 degrees of freedom. C: Linear regression for estimating the 
longer relaxation time 7^~ 1 . The t statistics assumes 445 degrees of freedom. 

4-1.2 Relaxation time for mean reversion 7 _1 



The relaxation time is also estimated by a linear regression over the logarithm 
of the empirical daily volatility autocorrelation function (11) given by: 

C(r) = ^kg£lHp fej±j _ i (23 ) 



We observe in figure 4B and 4C that IB3 data can be fit to two relaxation 
times: at r < 20 we fit 7f 1 = 20 ± 4 days and at r > 20 we fit 7 2 _1 = 270 ± 22 
days. 

The second relaxation time can be introduced into the original Heston model 
by making the long-term volatility fluctuate as a mean reverting process like: 



d9 = -7 2 [6(t) - 6 ] alt + K 2 ^e{t) dW 2 (t), (24) 

where dW 2 is an additional independent Wiener process. The autocorrelation 
function for the two time scales model acquires the following form [25]: 



"72T 

C(r) = -— + (25) 
cci a 2 

with a 2 = where 6 stands for the average of 6 given 6 . 
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The autocorrelation function fit in figure 4C indicates that relaxation times are 
of very different magnitudes and that it may be possible to solve the Fokker- 
Planck equation for such model approximately by adiabatic elimination. We 
pursue this direction further elsewhere. It is worth mentioning that another 
tractable alternative for introducing multiple time scales for the volatility is a 
superposition of OU processes as described in [26] . 

4-1.3 Amplitude of volatility fluctuations k 

From (5) the amplitude of volatility fluctuations is given by: 



The long term volatility 9 and the relaxation time 7 _1 have been estimated in 
the previous sections. As discussed in section 2.3, the constant a is related to 
the number of independent OU processes composing the stochastic volatility 
process as d — 2a. To avoid negative volatilities, a > 1 is required. We, 
therefore, assume a = 1 and calculate the amplitude of volatility fluctuations 
from (26) yielding « = 0.0032(10). 

4-1-4 Correlation between prices and volatilities p 

A nonzero correlation between prices and volatilities in (2) leads to an asym- 
metric probability density of returns described by (4). This asymmetry can be 
estimated either by directly computing the empirical distribution skewness or 
by calculating the empirical leverage function (13). Both procedures imply in 
computing highly noisy estimates for third and fourth order moments. In this 
section we propose estimating a confidence interval for p by computing the 
posterior probability p(p \ L), where L corresponds to a data set containing 
the empirically measured leverage function for a given number of time lags. 
Bayes theorem [27] gives the following posterior distribution: 



where Z(L) is a data dependent normalization constant. We assume igno- 
rance on the parameter p and fix its prior distribution to be uniform over the 
interval [— 1,+1], so that p(p) = U ([— 1,+1]). The maximum entropy prior 
distributions ^(7), p(9) and p(n) for the remaining parameters are gaussians, 
since their mean and variance have been previously estimated. The model 
likelihood for time lags < r < T < 7f 1 is approximately given by: 




(26) 




(27) 
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Fig. 5. A: Posterior density for the correlation between prices and volatilities p 
obtained by Monte Carlo sampling. We also plot in the same graph the cummulative 
distribution used to compute the 95% confidence interval indicated. B: Empirical 
leverage data as evaluated for the data set IB3. 



p(L | /o,7,0,k) = J 



dap(a) 

(27TCT 2 )* 



(28) 



x exp 



2o- 2 




M9- 



M 

£ 



X 



(1) 



OK 





where the first term inside the exponencial represents the empirical leverage 
function with being the M daily returns in the data set IB3. We choose 
p{o~) = U ([<7 min , o" max ]) to be uniform representing our level of ignorance on the 
acceptable dispersion of deviations between data and model. Having specified 
ignorance priors and the likelihood (28), we evaluate the posterior (27) by 
Monte Carlo sampling. In figure 5 we show the resulting posterior probability 
density and find the 95% confidence interval to be IC ^^y{p) = [—1.00, —0.42], 
what is strong evidence for an asymmetric probability density of returns. 



4-1.5 Probability density of returns 

Confidence intervals for the probability density of returns can be obtained 
by Monte Carlo sampling a sufficiently large number of parameter sets with 
appropriate distributions and numerically integrating (4) for each set. Dis- 
tribution features for each of the relevant parameters are summarized in the 
following table: 
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Fig. 6. Log- linear plot of empirical and theoretical probability densities of re- 
turns. Circles represent, from bottom to top, empirical densities at respectively 
1,5,20,40,80 and 160 days. Full lines indicate 95% confidence intervals obtained 
by numerical integration at Monte Carlo sampled parameter values. Densities at 
distinct time scales are multiplied by powers of 10 for clarity of presentation. 



Parameter 


Mean 


Standard Deviation 


e 


0.00101 days' 1 


0.00001 days' 1 




20 days 


4 days 




270 days 


22 days 


K 


0.0032 days' 1 


0.0010 days- 1 


P 


ICg 5% = [-1.00, -0.42] 


distributed as in fig. 5. 


a 


1.0 


fixed by theoretical arguments. 



We have independently sampled gaussian distributions on the parameters 
0,7i,/c and a uniform distribution U([— 1.00, —0.42]) for the parameter p. In 
figure 6 we compare empirical probability densities with theoretical confidence 
intervals at 95% finding a clear agreement at time scales ranging from 1 to 
160 days. 
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Fig. 7. Intraday autocorrelation function for each year composing the dataset IB2. 
The time scale separating micro and mesoeconomic phenomena is shown to be of 
about 20 minutes. 

4-2 High Frequency 

4-2.1 Autocorrelation of intraday returns 

Our main aim is to describe the fluctuation dynamics at intermediate time 
scales of formed prices (mesosconomic time scales) by a model which assumes 
uncorrelated returns. The price formation process occurs at time scales from 
seconds to a few minutes where the placing of new orders and double auction 
processes take place. We propose to fix the shortest mesoeconomic time scale 
to be the point where the intraday autocorrelation function vanishes. In figure 
7 we show that the intraday return autocorrelation function vanishes at about 
20 minutes for each of the 4 years composing data set IB2. We, therefore, 
consider as mesoeconomic time scales over 20 minutes. 



4-2.2 Effective duration of a day 

At first glance, it is not clear whether intraday and daily returns can be de- 
scribed by the same stochastic dynamics. Even less clear is whether aggrega- 
tion from intraday to daily returns can be described by the same parameters. 
To verify this latter possibility we have to transform units by determining the 
effective duration in minutes of a business day T e ff. This effective duration 
must include the daily trading time at the Sao Paulo Stock Exchange and the 
impact of daily and overnight gaps over the diffusion process. The Sao Paulo 
Stock exchange opens daily for electronic trading from 10 a.m. to 5 p.m. local 
time and from 5:45 p.m. to 7 p.m. for after-market trading, totalizing 8hl5min 
of trading daily. 

To estimate T e ff in minutes we observe that the daily return variance 

is the result of the aggregation of 20 minute returns, so that, considering a 

diffusive process, we would have: 



v 



(id) 



eff (20mm) 
20 



(29) 



15 




Fig. 8. A: Second cumulant of returns 02(f) versus the time lag in days for each 
year composing IB2 is used to estimate the mean daily variance v^ ld ^ following the 
procedure described in section 4.1.1. B: Plots of intraday second cumulants £2(5) of 
returns versus the time lag from 20 to 80 minutes employed to estimate, via linear 
regression, the 20 minutes variance v ( 20mm ) . Resulting T e ff 95% confidence intervals 
are also shown in the figure. 

It has been already observed that the volatility fluctuation dynamics are mean 
reverting with at least two time scales 7f 1 ^20 days and 72" 1 ~ 1 year. Con- 
sidering the longest relaxation time we estimate T e ff by estimating the mean 
daily volatility for each one of the years in IB2. In figure 8A we show lin- 
ear regressions employed for estimating the mean daily variance v ^ for each 
year in IB2 following the procedure described in section 4.1.1.. In figure 8B we 
show linear regressions employed to estimate y( 20mm ) ^ the effective duration 
confidence intervals are obtained from (29). The mean 95% confidence interval 
over the 4 years analysed results in IC^% (T e ff) = [9hl0min, 9h56min] what 
is consistent with 8hl5min of daily trading time plus an effective contribution 
of daily and overnight gaps. 



4-2.3 Probability density of intraday returns 

For evaluating the probability density of intraday returns we have reestimated 
the mean volatility 8 and the amplitude of volatility fluctuations k along the 
period 2001-2004 represented in the data set IB2. We then have rescaled the 
dimensional parameters as 9^ ID ^ = 6/T e ff, 7 (7Z) ) = ^/T e ff and k^ id ^ = n/T e ff. 
Having rescaled the distributions describing our ignorance on the appropriate 
returns we have employed Monte Carlo sampling to compute confidence inter- 
vals for the short time lags appoximation of the theoretical probability density 
described in (9). In figure 9 we compare the resulting confidence intervals and 
the data. We attain reasonably good fits for the longer time scales, as we ap- 
proach the microeconomic time scales the theoretical description of the tails 



16 



-0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 0.1 

X 

Fig. 9. Log- linear plot of empirical and theoretical probability densities of re- 
turns. Circles represent, from bottom to top, empirical densities at respectively 
20, 40, 80, 160 and 320 minutes. Full lines indicate 95% confidence intervals obtained 
by numerical integration at Monte Carlo sampled parameter sets rescaled by T e ff 
of the short t approximation for the probability density (9). Densities at distinct 
time scales are multiplied by powers of 10 for clarity. Note that the data shows tails 
fatter than predicted at shortest time scales. 

breaks with the empirical data showing fatter tails. 



5 Conclusions and Perspectives 

We have studied the Heston model with stochastic volatility and exponential 
tails as a model for the typical price fluctuations of the Brazilian Sao Paulo 
Stock Exchange Index (IBOVESPA). Prices have been corrected for inflation 
and a period spanning the last 15 years, characterized by memoryless returns, 
have been chosen for the analysis. We also have expunged from data a draw- 
down inconsistent with the supposition of independence made by the Heston 
model that took place in the transition between the first 22 years of long 
memory returns to the memoryless time series we have analysed. 

The long term mean volatility 9 has been estimated by observing the time 
scaling of the log-returns variance. The relaxation time for mean reversion 
7 _1 has been estimated by observing the autocorrelation function of the log 
returns variance. We have verified that a modified version of the Heston model 
with two very different relaxation times (77/ 1 ~ 20 days and 77 1 ~ 1 year) is 
required for describing the autocorrelation function correctly. We have used 
the minimum requirement for a non-vanishing volatility a > 1 to calculate 
the scale of the variance fluctuation k. Finally, we employed the Bayesian 
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statistics approach for estimating a confidence interval for the volatility-return 
correlation p, as it relies on a small data set to calculate a noisy estimate 
of higher order moments. The quality of the model is visually inspected by 
comparing the empirical probability density at time scales ranging from 1 day 
to 160 days, with confidence intervals obtained by a Monte Carlo simulation 
over the maximum ignorance parameter distributions. We have also shown 
that the probability density functions of log returns at intraday time scales 
can be described by the Heston model with the same parameters given that 
we introduce an effective duration for a business day that includes the effect 
of overnight jumps and that we consider the slow change of the volatility due 
to the longer relaxation time 

It is surprising and non-trivial that a single stochastic model may be capable 
of describing the statistical behavior of both developed and emerging markets 
at a wide range of time scales, despite the known instability and high sus- 
ceptibility to externalities of the latter. We believe that this robust statistical 
behavior may point towards simple basic mechanisms acting in the market 
microstructure. We regard as an interesting research direction to pursue the 
derivation of the Heston model as a limit case for a microeconomic model 
like the minority game [28]. For this matter we believe that the connection 
between the Heston model and Ehrenfest urns may be valuable. 

Perhaps the search for underlying symmetries and simple basic mechanisms 
that can explain empirical observations should be regarded as the main con- 
tribution of Physics to Economics. This contribution might be particularly 
useful to the field of Econometrics in which a common view is that a theory 
built from data 'should be evaluated in terms of the quality of the decisions 
that are made based on the theory' [29]. Clearly, these two approaches should 
not be considered as mutually exclusive. 
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